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There has been a new infusion of ideas in the study of the mechanism and early 
character of core-collapse supernovae. However, despite recent conceptual and 
computational progress, fundamental questions remain. In this all-too-brief con- 
tribution, I summarize some of the interesting insights achieved over the last few 
years. In the process, I highlight as-yet unsolved aspects of supernova theory that 
continue to make it a fascinating and frustrating pursuit. 



1 Introduction 



It has recently been shown that neutrino-driven Rayleigh-Taylor instabilities 
between the stalled shock wave and the neutrin^Dhcujs ( "Bcthe" convection) 
are generic feature of core-collapse supernovaeEiBBoOtj Whatever their role in 
reigniting the stalled explosion, their existence and persistence have altered the 
way modelers approach their craft. Supernovae must explode aspherically, and 
this broken symmetry is stamped on the ejecta and character of the blast, as 
well as on its signatures. Consequences of asphericity include significant grav- 
itational radiationPa natal kicks to nascent neutron stars0tl induced rotation, 
□ mixing of iron-peak and r-process nucleosynthetic products, the generation 
and/or rearrangement of pulsar magnetic fields, and, in extreme cases, jetting 
of the debris. 

However, there is no consensus yet on the centrality of overturn (or "con- 
vectionlllto the mechanism of th£_explosion itself, with some deeming it either 
pivotaipJa potentially importantpU'E3 or diversionary^ Nevertheless, all agree 
on the existence of convection in the gain region of the stalled protoneutron 
star, and this point must bfi-jstressed. A gain region is a prerequisite for the 
neutrino -driven mechanismt3 For heating to exceed cooling in steady-state ac- 
cretion, the entropy gradient must be negative, and, hence, unstable. Therefore, 
a gain region is always convective. In order to achieve quantitative agreement 
with the variety of observational constraints (explosion energy, residual neu- 
tron star masses, 56 - 57 Ni and tt N = 50" peak yields, halo star element ratios, 
neutron star proper motions, etc.), the "final" calculations must be done multi 
dimensionally. While if it can be shown that 1-D spherical models do explode 
after some delay, the true duration of that delay, the amount of fallback, and 
the energetics of the subsequent explosion must be influenced by the overturn- 
ing motions that can not be captured in 1-D. Convection changes not only 
the character of the hydrodynamics, but the entropies in the gain region and 
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the "efficiency:" ™. n f neutrino energy deposition that is the ultimate driver of 
the explosionlaOO Furthermore, implicit in a focus on 1-D calculations is the 
assumption that multi-D effects could only help, that they do not thwart ex- 
plosion. Hence, the belief that spherically-symmetric calculations are germane 
depends upon insights newly obtained from the multi-dimensional simulations. 
Nevertheless, it will be an important theoretical exercise to ascertain whether 
1-D models with the best physics and numerics can explode, if only because 
such has been a goal for decades. The "viablility" of 1-D models will be influ- 
enced in part by the transport algorithm employed (multi-group, flux-limited, 
full transport, diffusion), the microphysics (opacities and source terms at high 
and low densities), the effects of general relativity, the equation of state, psp^ 
vection in the inner core that can boost the driving neutrino luminositiesl^rr^l 10 ! 
and the inner density structure 



2 Neutrino Transport 

Though much of the recent excitement in supernova theory has concerned 
its multi-dimensional aspects, neutrino heating and transport are still central 
to the mechanism. The coupling between matter and radiation in the semi- 
transparent region between the stalled shock and the neutrinospheres deter- 
mines the viability and characteristics of the explosion. Unfortunately, this 
is the most problematic regime. Diffusion algorithms and/or flux-limiters do 
not adequately reproduce the effects of variations in the Eddington factors and 
the spectrum as the neutrinos decouple. Hence, a multi-group full transport 
scheme is desirable. 

To address the issues surrounding neutrino transport, we have recently 
created a neutrino transport code using the program Eddington developed 
by Eastman & PintoEj This code solves the full transport equation using the 
Feautrier approach, is multi-group, is good to order v/c, and does not employ 
flux limiters. The v e s, 9 e s, and V M s" are handled separately and coupling to 
matter is facilitated with accelerated lambda iteration (ALI). By default, we 
employ 40 energy groups from 1 MeV to either 100 MeV (D e and "i^s") or 230 
MeV (y e ) and from a few to 200 angular groups, depending on the number of 
tangent rays at the given radial zone, in the Feautrier manner. In this way, the 
neutrino angular distribution function and all the relevant angular moments 
(0'th through 3'rd) are calculated to high precision, for every energy group. 

The effect of the full Feautrier scheme vis-a-vis previouS'l!£E3QO'E!l cal- 
culations will soon be benchmarked and calibrated. However, we have already 
obtained several interesting results. Since the annihilation of v-v pairs into 
e + -e~ pairs depends upon the 0'th, l'st, and 2'nd angular moments of the neu- 
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trino angular distribution function, as well as upon the neutrino spectra, we 
can and have calculated the rate of energy deposition via this process exactly, 
though in the context of previous model runsu (still ignoring general relativ- 
ity). The v e + v e — > e + + e~ and v^ + v^ — > e + + e~ energy deposition rates in 
the shocked region are no more than 0.01 and 0.001, respectively, those of the 
dominant charged-current processes, v e + n — > e~ + p and v e + p — > e + + n. 
However, in the unshocked region ahead of the shock, depending upon the 
poorly-known ^-nucleus absorption rates, the v—v annihilation rate can be 
competitive, though it is still irrelevant to the supernova. These calculations 
should put to rest the notion that v—v annihilation is important in igniting the 
supernova explosion. 

It is thought that neutrino-electron scattering and inverse pair annihilation 
are the processes most responsible for the energy equilibration of the v^s and 
their emergent spectra. However, recent calculations imply that the inverse 
of nucleon-nucleon bremsstrahJung {e.g., n + n~^n + n + vv) is also impor- 
tant in equilibrating the z^'sEH This process has not heretofore been incorpo- 
rated in supernova simulations. Our preliminary estimates suggest that inverse 
bremsstrahlung softens the emergent v„ spectrum, since the bremsstrahlung 
source spectrum is softer than that of pair annihilation. In addition, given 
the large v^ scattering albedo, one must properly distinguish absorption from 
scattering, in ways not possible with a flux-limiter. Since the relevant inelastic 
neutral-current processes are stiff functions of neutrino energy, these transport 
issues beat directly upon the viability of neutrino nucleosynthesis (e.g., of n B 
and 19 F)B 

The new code allows us to calculate the difference bewteen the flux spec- 
trum (h v ) and the energy density spectrum (j„). The latter couples to matter 
and drives the supernova in the neutrino mechanism, while the former, or 
some variant of it, is frequently substituted for the latter in diffusion codes. 
Since matter-neutrino cross sections are higher for higher-energy neutrinos, 
the energy density spectrum is always harder than the flux spectrum. This 
hardness boosts the neutrino heating rates in the semi-transparent region. To 
illustrate this effect, in Figure 1 the ratio j v ^/h Ue is plotted versus neutrino 
energy at a time 30 milliseconds after bounce. The shock is then at 124 kilo- 
meters. It is clear that the ratio effect can be interesting. However, it is most 
pronounced in the cooling region below the gain region and tapers off as the 
shock is approached. Mezzacappa et aL0 in particular, have highlighted this 
correction, but self-consistent calculations from collapse to explosion, using 
the Feautrier or Boltzmann techniques (in principle equivalent), are needed, 
given the notorious feedbacks in the supernova problem. The same effect may 
be important in driving the protoneutron star wincP thought to be the site of 
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Figure 1: j v /h v versus for electron neutrinos 30 milliseconds after the bounce of a 15 M© 
core, using the code of Burrows & Pinto (1997) 

the r-process@0 Indeed, full transport calculations of r-process winds and 
the supernova, even in 1-D, will be illuminating. 

3 Conclusions 

In parallel with the ongoing evolution of the numerical tools being brought to 
bear on the supernova problem is the emerging realization that the systematics 
of the supernova phenomenon with progenitor is inching closer into view. As 
we unravel the mechanism, we simultaneously explore the origin of neutron 
stars and black holes, the birthplace of elements of which we are made, and 
the source of much of the energy of the ISM. As supernova modelers and the 
Jayhawks might say, ad astra per aspera. 
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